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Abstract 

We consider the effective two-Higgs-doublet potential with complex parameters, 
when the CP invariance is broken both explicitly and spontaneously. Diagonal 
mass term in the local minimum of the potential is constructed for the physical 
basis of Higgs fields, keeping explicitly the limiting case of CP-conservation, if the 
parameters are taken real. For special case of the two-doublet Higgs sector of the 
minimal supersymmetric model, when CP invariance is violated by the Higgs bosons 
! interaction with scalar quarks of the third generation, we calculate by means of the 

effective potential method the Higgs boson masses and evaluate the two-fermion 
Higgs boson decay widths and the widths of rare one-loop mediated decays H — ► 77, 
H - 99- 
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> : 1 Introduction 

{Sj . It is well-known that the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix originates 

from the Standard Model (SM) Lagrangian terms, describing the Higgs boson interaction 
with quarks (the Yukawa terms) 

O 

L = -gt^lHu^- gf^lHdi +h.c, (1) 

where = (u>, d') L , ffi = (&, s') L , = (t>, V) L , u£ = u' R , u% = c' R ,u% = t' R , d* = d' R , 
d\ = s' R , d R = b' R , and H denotes the scalar complex field doublet, Hk = euHl and gfj, 
gfj are the 3x3 matrices with matrix elements that are generally speaking complex and 
defined with an uncertainty coming from the phases of CP transformation 1 for the quark 
spinor fields and the Higgs boson scalar field. In order to diagonalize the quark mass 
term after spontaneous symmetry breaking H — > (0,v/y/2), the unitary transformations 
of the u % ' and d 1 ' quark fields u % LR = Ul,ru 1 [ r , d l L R = -Dl,_r d 1 ^ R are needed. After 
the diagonalization of the quark mass term the unitary matrices Ul and Di do not 
appear neither in the Yukawa Lagrangian terms (0) nor in the quark neutral current 
interactions, but arise in the quark u 1 ', d 1 ' charged current interaction terms gu' L , y^d' L W fl = 
guL r yfJJi J Didi J W tl . The product Vckm = UlD* l defines the complex CKM matrix, which 
decribes CP violation in the quark charged currents sector. In the framework of the 
SM the CP violation takes place since it is generally speaking not possible to get the 
mixing matrix with real matrix elements using CP transformations for six up- and down- 



-"^Let us remind, for example, that from the definition of the P transformation Pa^{p)P + = rj a a+ (—p) , 
where the complex factor \rj a \ =1 contains the P transformation phase, and a =0 or 1/2, it follows that 

P<p(x)P + =V 4>{x), Pip(x)P + =Vi/2 Joi/j(x'), where x = Px. 
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quarks. In other words, CP violation takes place in the SM because the number of quark 
generations is exactly three. 

There are other sources of CP violation besides the CKM mechanism. It is possible to 
introduce explicitly CP noninvariant hermitian Lagrangians [1 for the system of several 
scalar fields. For example, if we have three complex scalar fields ifi, <f2, <fi3 

L = XcpMrt + A V? W3 , CP L P + C + = L CP = Ae>^ 3 + X^tp^l 

where A is complex parameter and a is the CP transformation phase, not essential in this 
case. It can be rotated away by the phase transformation of the fields, related to charge 
conservation. One can see that L and L CP have different signs of the imaginary part 
of A. In this simple example the difference in the sign does not lead to any observable 
consequences, because the phase of A can be also rotated away by the U(1)q transforma- 
tion. However for the system with trilinear interactions of the four complex scalar fields 
it is generally speaking not possible to rotate away all phase factors. It is easy to show 
that the Lagrangian of such a system will be CP invariant only if the phases of the four 
parameters Aj respect certain conditions, which ensure the possibility to remove them by 
£7(1) rotations of the fields (pi. From this point of view the models with extended Higgs 
sector, where CP invariance of the Higgs potential with complex parameters is explicitly 
broken, are of particular interest. The simplest example is represented by the two-doublet 
Higgs potential of the MSSM, including (if the possibility of spontaneous CP violation 
[2] is not considered) ten parameters, four of them can be complex. In the framework of 
MSSM the dominant loop-mediated contributions from the third generation scalar quarks 
could lead to substantial violation of CP invariance of the two-doublet effective Higgs 
potential [3]. Various models with radiatively induced CP violation in the two-doublet 
Higgs sector have been studied [HE]. 

In this paper we develop further on our approach to the Higgs boson phenomenology 
in the scenario with CP violation considered in [S]. In Section 2, after brief introductory 
remarks, we calculate the effective Aj parameters of the two-doublet MSSM Higgs potential 
at the m top scale. In section 3 we consider in details the diagonalization of the mass 
term for the two-doublet Higgs potential with CP invariance broken both explicitly and 
spontaneously. In the Appendix some numerical results for the Higgs boson masses and 
the two-particle Higgs decay widths are presented. Our numerical results are compared 
with the output of other approaches. 



2 The effective two-doublet Higgs potential with CP 
violation 

In the general two-Higgs-doublet model (THDM) two SU(2) doublets of complex scalar 
fields are introduced: 

wm'^W'L r 2+ ) (3) 
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Their vacuum expectation values (VEV's) 



($1> = 72 ( I ) ' ( * 2) = 72 ( ) S 71 ( „«« ) ' (4) 

where v\ and v% are real. The phases (, relative phase of the VEV's, and £, relative 
phase of the SU(2) doublets, are introduced to consider the general case, their sum 9 
will be used for convenience of notations (section 3.3). For special case £ =0 the analysis 
of Yukawa sector with the two fermion generations can be found in [6] , where somewhat 
simpler form without the dimension 2 terms ${$ 2 + $ 2 $! an d real /z 2 2 , ^5,6,7 parameters 
of the THDM potential with spontaneous violation of CP invariance (( = 9 ^ 0) has 
been considered in the context of superweak (i.e. flavor- changing Higgs boson exchage 
mediated) CP violation in meson decays. 

The most general renormalizable hermitian SU(2) x U(l) invariant Lagrangian for the 
system of scalar fields (J2J), © can be written as 

C H = (V^V^t + (V U ^V^ 2 + K(V u ® i yV u $ 2 + k {V v ^V v ^x - U{$ 1} $ 2 ), 

(5) 

where 

U(* x , $ 2 ) = -/i?($I$ 1) - t4(& 2 $ 2 ) - ^2(^2)- /4 ($2$ 1) + (6) 
+A 1 ($ t 1 $ 1 ) 2 + A 2 ($5$ 2 ) 2 + A 3 ($ t 1 $i)($^ 2 ) + A 4 ($ t 1 $ 2 )($^i)+ 

+^($t$ 2 )($t$ 2 ) + ^($t$ 1 )($t$ l)+ 

+A 6 ($ t 1 $ 1 )($ t 1 $ 2 )+ A 6 ($I$i)($^i) + A 7 (^$ 2 )($ t 1 $ 2 )+ A 7 ($^ 2 )($ t 2 $i) 

The parameters fif 2 , A 5, A 6 A 7 are complex. Complex parameter k could be introduced 
to describe an interesting possibility of a mixing in the kinetic term [Zj. However, strong 
restrictions on the real part of k are imposed by precise experimental data on the gauge 
boson masses mw,z- Moreover, mixing in the kinetic term does not allow to construct 
the diagonal 4x4 matrix of the Higgs boson kinetic terms consistently with the diagonal 
matrix for their mass terms 2 . In the following we consider the case k — 0. 

Special case of the two-Higgs-doublet potential is the potential of the MSSM Higgs 
sector. At the energy scale Msusy (i-e. at the energy of the order of the sparticle masses) 
the tree level parameters A^...^ are real and can be expressed through the SU(2) x U(l) 
gauge couplings g\ and g 2 [E] 

M(Msusy) = HMsusy) = \{gl{M SUSY ) + gl(M SUSY )), (7) 
HMsusy) = \(gl(M S usY) - £ 2 (M 5C/sy )), K(M SUSY ) = - ±g 2 2 (M SUSY ), 
Xb(Msusy) = K(Msusy) = X 7 (M S usy) = 0. 

2 We analysed these conditions written in the form of ten linear equations, having the solution prac- 
tically only in the case n = 0. The mixed term is not obligatory to ensure the renormalizability. It is 
shown below that the contributions of self-energy diagrams absorbed by the Higgs boson wave-function 
renormalization to the effective parameters X^fij are zero, see also [H]. 
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At the scale Msusy the potential is CP invariant. However, the potential parameters 
of any model depend, generally speaking, on the energy scale where they are fixed or 
measured. The dependence is described by the renormalization group equations (RGE). 
The conditions ((ZJ) are the boundary conditions for the RGE. At the energies smaller than 
Msusy they are affected by large quantum corrections |H] where the main contribution 
is coming from the Higgs bosons - third generation quarks and scalar quarks interaction 
(the interactions with the first and second generations are suppressed). The potential of 
the Higgs bosons - scalar quarks interaction can be written in the form ^H] 



V° = V M + V r + V A + V 5 , (8) 



where 



V M = (-l^mj^ + Mj (QtQ) + M§U*U + M2D*D , (9) 

v r = rf ($}q) d + rf (i$f <t 2 q) u+ rf (gt$.) d*- rf (zQ V 2 $*) u* , (10) 



+ (11) 

A7j (®IQ) (Q'®j) + (^J°2® 3 )D*U + h.c] , z, j, k,l = 1,2, 

Vq denotes the four scalar quarks interaction terms, Pauli matrix er 2 = ( ^ . q ^ • The 

Yukawa coupligs for the third generation of scalar quarks are defined in a standard way 
h t = h b = ^2%. Following fTTl 3 : 

1 t) sin p ' u v cos p i ' 

rf 1;2} = % {-//*; At,}, r { ^ ;2} = ^ {A D ; -//*}, (12) 

they are complex in the case under consideration. One can observe CP violating terms 
of the structure similar to in the sector of Higgs-scalar quark interactions, so complex 
mixing matrices are expected to appear there. The trilinear parameters A t , A\> and the 
Higgs mass parameter /1 should be taken complex, the imaginary parts of the mixing 
matrix elements could be large. 

In the framework of the effective field theory approach ^U] the MSSM potential 
(JEJ) which explicitly describes sparticle interactions at the energy scale above Msusy is 
matched to an effective Standard Model-like Lagrangian at the energy scale below Msusy, 
where the sparticles decouple. So the MSSM effective Higgs potential at the energy scale 
m top , much smaller than M S usy, is represented by the general two-Higgs-doublet model 
potential (6), the parameters of the latter are expressed by means of the Higgs bosons - 
scalar quarks interaction parameters (|T2j) and the scalar quark masses, playing the role 
of ultraviolet Pauli- Villars regulators. The RGE boundary conditions (|ZJ) modified by the 
interactions of the third generation squarks with the Higgs bosons (these modifications 
are sometimes called the "threshold" effects, since the stops decouple at the M S usy scale), 
are imposed at the energy scale Msusy- They affect the evolution of Aj parameters, the 



3 For the case of CP conservation, considered in ^U], the trilinear parameters in (|10|) are real. Then 
r {i ; 2} = h u {-A*; A u}, 1^. 2} = h D {A D ; -fj] 
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Yukawa couplings h t) b and the gauge couplings g\^- We calculated radiative corrections 
to the boundary conditions (J7J) for \ parameters at the scale m top using the effective 
potential method [12]. The squark mass matrices (Aix) a b = 9 q V ^q* defined by (JEJ) were 
calculated and then substituted to the one-loop effective potential 



32tT S 



In 



a" 



decomposed in the inverse powers of Msusy- Taking into account the one- loop wave- 
function renormalization terms (i.e. terms introduced to absorb the contributions of 
self-energy diagrams to the Higgs bosons kinetic term, which are beyond the calculation 
by means of the effective potential method), the effective parameters can be evaluated as 
follows: 



Ai 



gl + gl 



+ 



32tt 5 



hi 



I A 



+ 



9i + gj 

4M S 2 USY 
A \{ ield + 



(h 2 t \(x\ 



M 2 



2 - 



I A 



6M 2 



SUSY 



h 



Ml 



6M S 4 USY 



+ 2hil + 



(13) 



+ 



768tt 2 



A 2 
A 3 



Ax (*<—>&), 



gl 



g\ 



4 
3 



1 - 



96n 2 M| USY 



16vr 2 



hi 



11^ + 9^-36 {gl + gl)ht)l 
{hi + hi) I 



I A 



+ 

12 



5 « 



/ + 2^ 



+ 



(14) 



A/f 2 

iw SUSY 



+ 



3(<7 2 2 -<7 2 ) [^(l/ip-IAI^ + ^d/xp 



I 3 " 

At?)] 



I A 



M 2 

JW SUSY 



+ 



128 vr 2 M s 2 USY 



+ A \{ ield 



Qgl 



ngf 



384 7r 2 



A, 



gl' 



+ 



I// 



167T 2 
2 



(A? + ^) ' 



/ + ^x tb 



96tt 2 M| USY L 



hi 



\A, 



M 2 

JKJ SUSY 



hi 



+ 



(15) 



M 2 

JW SUSY 



3fe 2 [^(|/i| 2 -|A| 2 ) + ^(H 2 -|A| 2 )] 



647r 2 M s 2 USY 



AA 



field 



where 



X, 



ft 



|A| 2 + |A| 2 + 2Re(AA) 
2-Msusy 



Mi 
64tt 2 ' 

- AAI 2 



M 2 

JKi SUSY 



6M S 4 USY 



The effective complex parameters A5 ; 6,7 
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Some details of the calculation can be found in [T^J. The one-loop wave-function renor- 
malization terms in (ll3l)-(fT5l) are 



A X f r ld = \{gl + g 2 2 )A' n , A \{ ie ld = \{g{ + g\)A 7 
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They are similar to the case of CP conservation containing the logarithmic contribu- 
tions and imaginary parameters as a consequence of ()12)) . and can be written as 
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Here and in the formulas given below I = In 



M, 



SUSY 



, where a = m top is the renormal- 



ization scale. The one-loop wave-function renormalization does not yield a CP violating 
contribution to Aj. For convenience we introduce the notation for the deviations of effec- 
tive parameters Aj from Xf USY = X^Msusy) following [3]: 



Ai,2 = Ai^C — AAi^/2, A3 t 4 = A 
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finite 



A A{^7 = ■ 



(22) 

(23) 
(24) 



In the end of this section we would like to make some general comments as well 
as some comments in connection with results obtained by other authors. Like in the 
existing effective field theory approach JU] we are using the standard scheme of leading 
logarithmic terms resummation by means of RGE, additionally taking into account in the 
boundary conditions at the scale Msusy the effects of Higgs bosons - third generation 
of scalar quarks interaction. The one-loop effective parameters (j!3)l - (fT9*|) satisfy the 
boundary conditions defined by and modified by the soft supersymmetry breaking 
potential terms (jHJ) ("threshold effects"). The terms with the logarithmic factor I describe 
the parameters evolution from the energy scale Msusy down to the scale a = m top . 
Finite power term threshold corrections to Ai 7 appear from the so-called F-terms (the 
trilinear interaction terms in (fTUj) ) and D-terms (contained in (fTTf ). The corrections 
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l-loop+2-loop + l-loop(D+wfr) 


0.715 
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-0.360 


0.033 



Table 1: Numerical comparison of various corrections to the A^ parameters at the scale 
mtop- For convenience of the following Higgs boson masses comparison, the same pa- 
rameter values as in the package CPsuperH are chosen here: rriz = 91.19 GeV, 
m b = 3 GeV, m t = 175 GeV, m w = 79.96 GeV, g 2 = 0.6517, g x = 0.3573, v = 245.4 GeV, 
Gp = 1.174- 10~ 5 GeV- 2 , a s {m t ) = 0.1072, tan/? = 5, M SUSY = 500GeV, a = m t , 
m H± = 300, \A t \ = \A b \ = A = 1000 GeV, = 2000 GeV, <p = arg(/iA ii6 ) = 0. The 
abbreviation 'wfr' stands for the 'wave-function renormalization'. 



to A5 come from the F-terms only. Radiative corrections to the parameters Ai...j of 
the effective two-Higgs-doublet potential have been considered earlier in [3] for the case of 
broken CP invariance and in fl] for the case of CP conservation. Phenomenological 
consequences of the two-doublet system are usually analysed assuming for simplicity A t = 
A b and introducing the universal phase axgnA tbl so that A 5 = |A 5 | exp [i 2 arg(/x/l)], 
A 6 = I A 6 | exp [i arg(/M)], A 7 = | A 7 | exp [i arg(/M)]. 

Only the leading D-term contributions were calculated in [Hj, [0|. m our expressions 
for the effective parameters (fT3*|) - (fTn|) the nonleading D-term contributions are represented 
by the power terms containing gauge couplings gf, g\. The one- loop contributions of the 
wave-function renormalization AA{ iei< 4 are neglected in [3], [Tl]. However, the QCD and 
weak corrections to Yukawa couplings up to two loops, not calculated in our case, have 
been included there. The expressions for A 1,2,3,4 (fH^ - ffTKj) do not contain imaginary parts 
up to the two-loop approximation and coincide with the results of [Sj, [Ej if we omit the 
contributions of nonleading D-terms and AA{"^ \ terms. If \i and A are real, the expres- 
sions (fl"3*j) -(fl"9" j) are consistent with the results of [TU], where the .D-terms contribution was 
calculated 4 . Let us note that it is not possible to generalize the expressions for real A 5i6i 7 
in the case of CP violating potential by the straightforward replacement of the real fi, A 
parameters to the complex ones. 

If we neglect the contributions of .D-terms, the wave-function renormalization terms 
AA{ 2e? 4 and the terms of the order of h% for the 6-quark couplings, only the one-loop 
corrections of the order of O(hf) remain. This approximation was discussed in [T4*l I12j. 
For example, A2 is given by 



g 2 2 + g\ 



32 7T 2 



\Af ( \A\ 2 \ { 



(25) 



1 In (|13(l - (|15|l we kept the terms of the order of gf 2 - 
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The beta-function for A2 contains large negative contribution — 6/tf [H]\, or equivalently, 
A2 (13) contains the large logarithmic term 6hfl /(32it 2 ) which was observed in the first 
calculations [9 . In the following the negative AA 2 defined by (22) gives large positive 
contribution to the light Higgs boson mass in (38). 

Numerical comparison of the A; parameters evaluated using different approximations 
is presented in the Table 1, where for our case in the second line of the Table 

A A i = {one-loop contribution} + {one-loop(P — terms + wave-func.renormalization)}. 

One can conclude that the one-loop corrections from .D-terms and wave-function renor- 
malization can be of the order of the leading two-loop corrections. Difference of the 
effective A» of the order of 10 _1 may result in the deviation of Higgs boson masses around 
5 GeV and even more. 



3 Diagonalization of the effective potential mass term 
in the local minimum 

3.1 Complex n\ 2 , ^5,6,7 parameters, 9 =0 

The components Ui, r/i, Xi of the SU{2) doublets ©, © are not a physical Higgs fields 
(mass eigenstates). In order to extract the Higgs boson masses and the self- interaction 
of the physical fields from the potential (J2J) it is necessary to diagonalize the mass term 
of the latter in the local minimum. This problem has been considered in |H] for the 
case of complex n\ 2 , As^j parameters and the zero phase of the $2 VEV 9 =0. The 
diagonalization of the mass term is performed in two stages. First the CP-even fields 
h,H, the CP-odd field A ('pseudoscalar') 5 and the Goldstone field G° are defined by the 
linear transformation 

h = — 771 sin a + 772 cos a, (26) 
H = r/i cos a + i]2 sin a, (27) 
A = -xisin/3 + X2Cos/3, (28) 
G° = xicos/5 + x 2 sin/3, (29) 

where tg/5 = v 2 /vi and 

S2fs(m 2 A + m 2 z ) + tj 2 ((AA 3 + AA 4 )s 2/ 3 + 2c%ReAA 6 + 2s%ReAA 7 ) 
S a ~C2f}{m\ - m 2 z ) + v 2 (AX 1 c 2 p - AA 2 s| - ReAA 5 c 2/ 3 + (ReAA 6 - ReA\ 7 )s 2 p) ' 

Here the relations g\ + g\ = g 2 m 2 z /m 2 v , g\ — g\ = g\ (2 — m|/m^) are used. Then we 
substitute to the effective potential the real parameters Ai.,2,3,4 and the real parts 
ReyU^ 2 , ReAs^,?, which are related by linear transformation 0113 HE]: 

Ai = ^i(-) 2 m 2 h +(^) 2 m 2 H -^Rerf 2 ] + i(ReA 7 tg 3 /3 - 3ReA 6 tg/3), (31) 

2V Z C/3 C/3 &p 4 



3 The fields h,H,A are the physical fields at <p = aig(fiAt t b) = 0, rm. 
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A 2 = -^[(^) 2 m2 + (^) 2 m^-^Re / iy + i(ReA 6 ctg 3 /3-3ReA 7 ctg/3), (32) 
2v 2 sp sp sp 4 



\ 1 r 2 Re A i 12 , S 2a, 2 2m Re ^6 . a ReA 7 

A 3 = — [2m H± + — (m H - m h )\ — ctg/3 — tg/3, (33) 

spcp s 2 p 2 2 



1 ,Re/i 2 2 2 2 ReA 6 ReA 



A 4 = ^(^^ + m 2 4 -2m 2 ^ ± )-=ctg/3-^tg/3, (34) 



1 ,Re/ij 2 2\ ReA6 „ ReA 



ReA 5 = ^(^H- m ^)-^p c tg/5-^tg/?, (35) 

f S/?C^ 2 2 

/i 2 = Aiw 2 + (A3 + A4 + ReA 5 )^-Re/i 2 2 tg/5+^(3ReA 6 ctg/3 + ReA 7 tg/5),(36) 

l4 = A 2 ^ 2 2 + (A 3 + A 4 + ReA 5 )y - Re^ 2 2 ctg/3 + -^(ReA 6 ctg/3 + 3ReA 7 tg/3).(37) 

At the purely real parameters (in the following we shall name this case of ip = as the 
CP-conserving limit, ReAj = |Aj|, ReAAj = |AA;|) the relations flSEJ), (}3T|) set to zero the 
potential terms which are linear in the fields, so they are the minimization conditions. It 
follows from the equations (|31 )) -(|35 j) that in the CP conserving limit the CP-even Higgs 
boson masses and the real part of the fi 2 2 parameter can be expressed as 

m l = s 2 +/3 m| + c 2 a _ p m 2 A - (38) 
-w 2 (AA lS 2 c 2 + AA 2 c 2 s 2 - 2(AA 3 + AX A )c a cps a sp + ReAA 5 (s 2 s 2 + c 2 cj) - 
-2c Q+/3 (ReAA 6 s a c /3 - ReAA 7 c Q s /3 )), 



m 



2 2 2,2 2 tor\\ 

H = c a+/3 m z + s a _pm A - (39) 
-v 2 (AA lC 2 4 + AA 2 s 2 4 + 2(AA 3 + ^)c a c p s a s p + ReAA 5 (c 2 4 + s 2 a c 2 p) + 
+2s a+/ g(ReAA 6 c Q c / 3 + ReAX 7 s a Sp)), 
v 2 

m H± = m w + m A~^ (ReAA 5 -AA 4 ), (40) 

v 2 

Re/x 2 2 = spCp[m A - — (2ReAA 5 + ReAA 6 ctg/9 + ReAA 7 tg/3)]. 

After the substitution of (|3*T|) - (j33j) . (jHSI), fj3*7j) to (0) we find the mass term of the effective 
potential 

U mass (h,H,A) = coA + Cl hA + c 2 HA+ (41) 

2 2 2 H 

The minimization condition Co =0 fixes the imaginary part of the /z 2 2 parameter 

2 

and the factors in front of the nondiagonal terms hA and HA in the local minimum Co = 
have the form 

v 2 

ci = y(s Q S/3 - CaC^ImAs + v 2 (s a cplm\ 6 - c a s /3 ImA 7 ), (43) 
v 2 

c 2 = --^-(sqC/j + c a sp)lm\ 5 - v 2 (c a cplm\ 6 + s a s p ImX 7 ). 



Im/i 2 2 = — (s/jc^ImAs + ImA 6 + ImA 7 ), (42) 



They include only the imaginary parts of the parameters Im/i^, ImAs^j. The nondiagonal 
term hH does not appear in (j41j) . so in the mixing matrix ()45|) M 12 = M 21 = 0. 

At the second stage in order to remove the nondiagonal terms hA and HA we perform 
the orthogonal transformation in the h, H, A sector 





(h, H, A) M \ H \ = (h u h 2 M)ai k M 2 kiaij \ h 2 \ , (44) 
where the mass matrix is 

l ( m 2 h ci 
M 2 = - m\ c 2 | (45) 
\ ci c 2 

and get the physical Higgs bosons hi, h 2 , h% without a definite CP parity 6 . The eigen- 
values of the M 2 matrix define their masses squared and the components of normalized 
eigenvectors are the matrix elements in the rows of the mixing matrix a^. The squared 
masses of Higgs bosons are {m 2 hl < m\ 2 < m\ 3 ) 

6 + 27t\ a 2 

3 ) ~ ~3 

6 + 47r\ a 2 



< = 2 VHz) cos ( — - ^ , (46) 



m l 2 = 2 \/(-v) cos 



where 



m h, = 2 v /(-g)cos 



9 = arccos 



3/3 
3j 3 ' 



-q 



31 



r = F7( 9a i fl 2 - 27a - 2a\) , q = J-(3ai - a\) , 
54 9 

22,22,2222 2 2 2 

ai = m h m H + m h m A + m H m A — c l — c 2 , a 2 = — m h — m H — m A , 

2 2,2 2 2 2 2 

a = c 1 m H + c 2 m ft — m h m H m A . 
The normalized eigenvector components (h,H,A) = a^hj, = a'ij/nj are given by 



%i = (( m H ~ r n hl )(m A - m hl ) - c 2 ), a 21 = c x c 2 , a 31 = -c^m^ 



a 12 = -cic 2 , a 22 = -{{ml - m 2 h2 ){m 2 A - m 2 h2 ) - c\), a' 32 = c 2 {m\ - m 2 h2 ), 

' l 2 2 \ ' I 2 2 \ ' I 2 2 \l 2 2 \ 

13 = -ci(m H - m h3 ), a 23 = -c 2 {m h - m h3 ), a 33 = {m h - m h3 ){m H - m h3 ), 



rii — y {a 2 i + a 2i + a 3i ). The Higgs boson masses m^, rrih 2 , mh 3 and the mixing matrix 
elements a^, which describe the mixed states, are shown in Fig.2-4 as a function of the 

6 Note that this picture is different from the well-known description of weak CP violation in meson 
decays, when the mass splitting Am of the states is given by 2ReAii2, Mn the off-diagonal elements of 
the complex 2x2 mass matrix, and the meson mixing e parameter is ImMi2 / (V2Am) . The meson decay 
formalism uses the non-hermitian effective Hamiltonian and not precisely orthogonal mass 'eigenstates'. 
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A ti b, n parameters and/or the universal phase <f = arg(/nA t> {,). Different to the figures 
in [5J, the m#±, tg/3 parametrization is used for the convenience of comparison with 
and The parameters c\ and c 2 can change a sign with the variation of the phase (p, 
the ranges of positively or negatively defined c\ and C2 depend on the primary choice of 
the m H ±, tg/3, A, /i and Msusy in the CP conserving limit. When we pass the zeroes 
of c\ and c 2 , the matrix elements ay are expected to change their signs respecting the 
requirement of the left orthonormal basis for the eigenvectors. It is essential that m^, 
m/j 2 and are positioned in the mass matrix along the diagonal from the upper left 
to the lower right corner, satisfying in the limiting case c\ = c 2 = the correspondences 
— > min(mh,mH,mA), rn^ — > max{my l ,mH-,' m A) ("the mass ordering"). Note also 
that as AAj increases, the denominator of (jHOj) can change sign, so for the mass ordering 
one must define the angle oi(cp) consistently with the boundary condition at the scale 
Msusy, which has the known form m A + m 2 z = — sin2a/sin2/5(m^ — m\), following from 
dSU-dMj) and ©. 

Some numerical values for the Higgs boson masses , rrih 2 , rrih 3 as a function of the 
phase if in our approach, and masses of the states Hi, H 2 , H 3 evaluated by means of 
CPsuperH [I7| and FeynHiggs [18J packages are shown in the Table 2. These packages 
are using the renormalization group improved diagrammatic calculation that uncludes 
radiative corrections to Yukawa couplings up to two-loops. Detailed general discussion 
on the conciliation of results obtained in the frameworks of the diagrammatic and the 
effective field theory approaches can be found in ^Hj- Different renormalization schemes 
in which calculations in the two approaches are performed, may lead to the deviations of 
results evaluated with parameters taken at different renormalization scales, so the untrivial 
reevaluation of parameters is needed for consistency. Besides this it is important to notice 
that in the CPsuperH and FeynHiggs packages the SU(2) eigenstates 771,2, £1,2 are directly 
transformed to the Higgs boson mass eigenstates, which is different from our procedure, 
when we first transform to the states of the CP-conserving limit and then rotate to ^1,2,3- 
The 'intermediate' Higgs boson states (h, H, A) of the CP conserving limit are not used, 
so the 771, r]2 mixing angle a is not introduced there. For this reason at if =0 the analogue 
of the mixing matrix (44), has nonzero off-diagonal matrix elements ai2 = a 2 i 7^0 

and in the analogue of the mass matrix (47) mu and 17121 (the hH mixing terms in our 
notation) are also nonzero. In the framework of the 'direct' diagonalization procedure the 
matrix elements of (45) have the form 

rriu = % Sp + v Re\ 5 + v Re\ 6 s 2 /3 + 2v XiCp, 
m 2 2 = m 2 A c 2 p + v 2 Re\ 5 c 2 p + v 2 Re\ 7 s 2l 3 + 2v 2 \ 2 s 2 p, 
777,12 = f 2 ReA 6 c 2 p + sp (v 2 Re\ 7 S[3 + cp {—m 2 A + v 2 A 3 + f 2 A 4 )), 

mi3 = -iw 2 (2 ImA 6 C/3 + ImAsS/3), 

m 23 = --v 2 (lm\ 5 C/3 + 2 ImXjSp), 
"7.33 = rn 2 A 

and the parameters do, ai, a>2 in (46) should be redefined as follows 

a = m 2 2 m 33 + 77723 m n + m 2 3 m 22 - 2 mi2 m 2 3mi3 - m n m 22 m 33 , 
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Table 2: The Higgs boson masses (GeV) in our case and calculated by the packages 
CPsuperH [iTj and FeynHiggs 18] (in the one-loop regime) at the same parameter val- 
ues a EM (m z ) = 0.7812-10" 2 , a s (m z ) = 0.1172, G F = 1.174 ■ 10~ 5 GeV -2 , tan/5 = 5, 
Msusy = 500 GeV, \A t \ = \A b \ = A, = 2000 GeV, A = 1000 GeV, m H ± = 300 GeV. 

at = m u m 22 + m u m 33 + m 22 m 33 - mf 2 - mf 3 - ml 3 , 
a 2 = - m u - m 22 - m 33 

We checked that both the 'two-step' and the 'direct' diagonalization methods lead within 
our procedure, as expected, to the same masses of Higgs states m^, m^ 2 , m/j 3 (see Table 
2). For the parameter values in the comparison, Table 2, the benchmark point of the 
maximal CP violation 'CPX scenario' 20] at Msusy =500 GeV was used. Extended list 
of numbers (Table 5) including also the rare one-loop mediated decay widths hi — > 77, 
hi — > gg and the tree-level two-particle decays hi — > ff can be found in the Appendix. 
Good qualitative agreement of results is observed, but diversity of approaches to the 
calculation of radiative corrections makes precise numerical comparisons difficult. 

3.2 Real /i^ 2? A5 6,7 parameters, 6 ^0 

If the parameters /x 2 2 , As^j of the effective potential (j2j) are real, the latter is CP invariant. 
It is easy to show 013112], that the phases of complex parameters /i 2 2 , As^j can be rotated 
away by the U(1)y hypercharge transformation if the conditions 

Im(/4 A* 5 ) = 0, Im(/i 2 2 A 6 ) = 0, Im(/i 2 2 A 7 ) = 0. (47) 

are satisfied. Insofar as the physical motivation of these 'fine tuning' conditions is not 
available, the case of real parameters and nonzero phase 9 of the VEV, when CP is broken 
spontaneously, looks rather artificial. The local minimum of the effective potential (0) 
occurs at A 5 > (i.e. purely imaginary fiA, see (|I7j) ) and 

2 vf •> v\ \ 

cosg= ^ 12 ~~ T . (48) 

\ h ViV 2 
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Combining this equation with the diagonalization condition (jHHjl we get 

2 

cos0=^ + l, (49) 
A 5 v 2 

so there is no minimum if m 2 A > 0. In the case A 5 < (|48j) corresponds to the maximum, 
the absolute minimum is achieved at the endpoints cos# = ±1. For example, the absolute 
minimum at 9 =0 (taking into account again the diagonalization condition (jSSJ)) is absent 
if 

m 2 A >2\\ b \v 2 . (50) 

and it follows that for the case of real /z 2 2 , A5 6J and CP broken spontaneously there are 
no mass eigenstates in the framework of our diagonalization procedure, at least if rriA is 
not extremely small. 



3.3 Complex /if 2 ? A5 6,7 parameters, 9 ^0 

In the case of complex parameters and the nonzero phase of $2 vacuum expectation value 
7 , the CP invariance of the potential is broken both explicitly and spontaneously. The 
condition to set to zero the derivative dU/89 includes both the real and the imaginary 
parts of fj, 2 2 and A 5i6i7 : 

cos #(2Im/i 2 2 — v 2 ImA6 — v \ ImA 7 ) — Uii^ImA 5 cos 29 + (51) 

+ sin 9{2Reji\ 2 — v 2 ReA@ — v 2 Re\ 7 ) — ViV 2 Re\ 5 sin 26 = 0. 

The condition of the extremum for Im/i 2 2 depends on the phase between the VEV's 9, 
while the diagonalization condition for Re/i 2 2 depends also on the relative phase £ (see 
©, ©) of the SU(2) doublets. At the real //f 2 , A 5i6i7 and 9 ^ the equation (JHTJ is 
reduced to lJ35)l . 

For convenience we present the extremum conditions dU/dr] = 0, dll/dC, = in the 
cases of zero and nonzero 9 in the form of Tables 3 and 4, where the factors in front of 
the potential parameters are shown. Bulky condition for the real part of fx\ 2 to define the 
pseudoscalar mass itla for the general case of nonsero phases can explicitly be evaluated 
as follows: 

Re/4 = (52) 

_ A v 2 cos 9 sin 3 (2(3) sin 2 (fl + Q) 

2 3 + ( 1 - cos 9 cos (cos 4 p - f sin 2 (2/3) ) + sin 4 (3 + cos 9 cos £ ( 1 - sin 4 0) + 



ImA 1 



w 2 (cos 4 (3 cos 2 £ + cos 2 9 sin 4 (3 + cos (3 cos(9 — £) sin/? sin(2/3) 
cos 2 (3ctg(3 sec + cos £ sin(2/3) + sec sin 2 /3tg/3 

u 2 (sin 2 (2/3) sm(9 - Q + sin 4 /?(sin(2fl) + tgg) + cos 4 /3(tgfl - sin(2Q) 
2(cos 2 /3ctg/3 sec + cos £ sin(2/3) + sec 6 1 sin 2 (3tg(3) 

+ ReA 6 i^ 2 cos 2 (3 + 



7 The upper component of ($2) in Q is taken to be zero. Otherwise additional constraint for the VEV 
components should be imposed to ensure the existence of the massless gauge field (photon) [23 
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9 = 


Ai 


v\ 










A 2 








v 2 




A 3 


2 U 2 


2^2 


¥1 




A 4 


UP- 


2 W 2 




2^1 


ReA 5 


I,.2 
2 U 2 


2 U 2 


2 U 1 


2 U 1 


ImAs 


-\v\\g9 





-iu?tg0 





ReAg 


\v x v 2 {2 + cos 29) sec 6 


\v x v 2 


|t) 2 sec 9ctgj3 


|w 2 ctg/? 


ImA6 


—V1V2 sin 6 1 











ReA 7 


|f 2 sec 0tg/3 




\v x v 2 {2 + cos 26>) sec6» 


\v x v 2 


ImAy 








— f iv 2 sin # 





Refx 2 2 


— tg/3 sec 9 




— ctg/3 sec 9 


-ctg/3 



Table 3: The factors of the extremum conditions for /i 2 fi 2 at zero and nonzero 9. 

v 2 cos 3 /3 sin (3 sin £ 



+ ImA 
+ReA 7 I 



+ 



+ImA' 



' cos 2 (3ctgf3 sec 9 + cos £ sin(2/3) + sec 9 sin 2 /3tg/3 

w 2 cos 4 /3(4 cos(0 + 2^-2 cos(2#) sec fljtgff 
4(cos 2 /3ctg/3 sec 6> + cos £ sin(2/3) + sec 9 sin 2 /?tg/3) 

v 2 (2 sin 2 (2/3) cos£ + 2 sec sin 4 /? - cos(2# + £) sin 2 (2/?))tg/f 
4 (cos 2 /3ctg/3 sec 9 + cos £ sin (2/5) + sec # sin 2 /3tg/5) 

v 2 sin(2/?)(2 cos 2 /3 cos f sin (0 + + sin 2 /3(2 sin£ + sin(2fl + 0)) 
2(cos 2 (3ctg(3 sec + cos f sin(2/3) + sec # sin 2 /3tg/?) 

_ Im// 2 sin(2/3) sin£ + 

12 cos 2 PctgP sec 6 1 + cos £ sin(2/3) + sec 6 1 sin 2 /3tg/3 

^ 2 \ 



cos 2 f3ctg(3 sec + cos £ sin(2/3) + sec 6* sin (3tg(3 

If we set 9 = and £ = 0, the formulas coincide with the special case of only the explicit 
CP violation (|35|h (}4*2*|) . The substitution of the extremum conditions corresponding to 
Tables 3 and 4 to (J5T|) gives an identity independently on the expression (|52~|l for Re/i 2 2 . 
The extremum is a minimum if the second derivative in 9 is positively defined 

— sin 6 l (2Ini / u 2 2 — t> 2 ImA 6 — t^ImAy) + 2w 1 t> 2 ImA 5 sin 29 + (53) 
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Table 4: The factors of the extremum condition for Re/i 2 2 at 9 = and for Im/i 2 2 for zero 
and nonzero 9. 

+ cos 6'(2Re / u^ 2 — v 2 Re\ 6 — t> 2 ReA7) — 2viV 2 Re\ 5 cos 29 > 0. 

Numerical investigation shows that this condition is fullfilled in a rather wide range of 
the MSSM parameter space. If for simplicity we set £ =0 then the second derivative is 
positively defined in any region of the parameter space, so no restrictions on the phase of 
spontaneous CP breaking appear in this special case from the minimization. 

The diagonalization of the effective potential mass term in the local minimum for 
the general case 9 ^ and £ ^ is performed analogously to the procedure described 
in section 3.1 using the following scheme: (1) we define the four h, H, A, G° linear 
combinations of independent fields rji, rj 2 , Xu X2 that are contained in the two-doublet 
system (J2j), (JHJ), where for the Goldstone field G° we define zero row of matrix elements 
and zero column of matrix elements in the symmetric mass matrix 4x4. In other words, 
the Goldstone mode is introduced as the linear combination, orthogonal to the plane 
defined by the "directions" in the complex scalar fields space, parallel to the VEV's V\ 
and v 2 exp + ()}. Then the mass matrix 4x4 includes the symmetric 3x3 block with 
zero matrix elements in the power of the extremum conditions from Tables 3 and 4; (2) 
we define an orthogonal transformation for the 3x3 submatrix fixing the mixing angle a 
in the sector h — H to set to zero the hH nondiagonal term. In the framework of this 
procedure for the case of nonzero phases £ 7^0, 9^0 (when the fields are denoted by the 
symbol " ~" ) the limiting cases of zero phases £ = 9 =0 (when the notation for the fields 
does not contain the symbol " ~") and also the CP conserving limit in the mass basis h, 
H, A, are clearly seen. For the physical Higgs fields in the case £ = 0, 9 7^ we finally 
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obtain the representation 



h = — rji sin 5 + (x2 sin 6* + 772 cos6*) cos 5, (54) 

H = rji cos 5 + (%2 sin# + 772 cos^) sin 5, 

A = — Xi sin /3 + (X2 cos 9 — rj 2 sin 0) cos/3, 

G° = xi cos/3 + (x2 cos 9 — X2 sin 9) sin/3. 

We checked explicitly, using the symbolic calculation packages, that direct substitution 
of these fields to the potential (J2J) gives the symmetric 4x4 squared mass matrix with 
zero row and column, corresponding to the Goldstone mode. The non-diagonal matrix 
elements of the 3x3 block, corresponding to the nondiagonal terms hA HA in the local 
minimum, can be written in the form 

v 2 

c\ = — — (cos(a+/3) cos(2^)Im A5 — 2 sina cos/3 cos^Im A6 + 2 cosa sin/Jcos^Im A7— (55) 



cos(a + (3) sin(26>)Re A 5 — 2 sin a cos /3 sin 9Re A 6 + 2 cos a sin (3 sin 9Re A 



7), 



v 2 



C2 = — — (sin(a + /3) cos(2^)Im A5 — 2 cosa cos/? cos^Im A6 + 2 sina sin/3 cos^Im A7+ (56) 

+ cos(<5 + (3) sin(26 l )Re A 5 — 2 cos a cos (3 sin^Re A 6 + 2 sin a sin (3 sin 9Re A 7 ). 

In the case 9 = they coincide with (|43|) . 

The same scheme is suitable for the case £ 7^0, 9 7^0 when the relative phase £ be- 
tween the SU{2) doublets appears in the mass eigenstates, which are obtained by the 
replacement 9^9 — ^: 

h = — rji sin a + (X2 sin(6 l — £) + 77 2 cos(9 — £)) cos a, (57) 

H = 771 cos 5 + (%2 sin(6> — £) + 77 2 cos(6* — £)) sin 5, 

A = -Xisin/3+ (x2cos(0 — £) — r) 2 sxn(0 — £)) cos/3, 

G° = xiCos^+ (x 2 cos(6>-£) -X2sin(6»-£))sin/3. 



4 Summary 

The potential of a two-Higgs-doublet model in the general case is not CP invariant and 
the parameters fi 2 2 and A 5;6)7 of the two-doublet MSSM Higgs sector should be taken 
complex. The choice of purely real parameters implicitly assumes that the fine-tuning 
conditions (|4*7|) are additionally imposed without clear physical motivation. In the MSSM 
the complex parameters naturally appear if we allow the CP invariance violating mixings 
in the squark-Higgs boson sector of the MSSM, analogous to the CKM mixings for the 
three quark generations in the charged current sector of the Standard Model. If these 
mixings lead to a strong CP parity violation 8 and the scalar sector of the MSSM is cou- 
pled strongly enough (i.e. large imaginary parts of the parameters /jf 2 and ^5,6,7 appear), 
the deviations of the observable effects in the scenario with CP violation from the phe- 
nomenology of the standard scenario can be substantial. The deviations are particularly 

8 Recent discussion of the weak CP violation scenarios can be found in |22j . 
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strong if the power terms A^/Msusy, ^/Msusy are large and the charged Higgs boson 
mass does not exceed 150-200 GeV, being rahter weakly dependent on the value of tg/3. 
Such models could lead in principle to a reconsideration of the experimental priorities [2H] 
for the signals of Higgs bosons production in the channels 77, bb, W + W~ , ZZ, ttH, bbH 
etc. at the LHC. The scenario with light Higgs boson ~ 70 — 80 GeV that could 
escape the detection at LEP2 [21], the analysis of h± signal at Tevatron and the high- 
luminosity linear colliders [23] demonstrate that physical possibilities in the framework of 
CP violating scenarios could be considerably modified in comparison with the traditional 
CP conserving limit. 

The comparison of our results for the masses of scalars m^, rrih 2 and and their 
two-particle decay widths with outputs of the CPsuperH [T2j and the FeynHiggs [TH] 
packages demonstrates rather good qualitative agreement. However, is some cases high 
sensitivity of the observables to the magnitude of radiatively induced correction terms in 
the effective two-Higgs-doublet potential shows up, so careful complementary analysis of 
the theoretical uncertainties is appropriate. 

The relative phase of the SU{2) scalar doublet ( and the VEV phase £ (J1J) could 
be constrained on the basis of the conditions for the mass term diagonalization and the 
potential minimization (Section 3.3). In principle these conditions could lead to some 
nontrivial relations between the (, £ and the variables of the MSSM parameter space. 
However, at the first sight it is questionable to expect some direct relations of this type 
connecting the CKM phase and the (, £ phases of the THDM, which seem to describe 
the CP violation of different origin. Returning to the notations of the Introduction, we 
can write the THDM type II Yukawa term as 

— L = r ? «^$ 1 + ^.^4$ 2 + h.c., (58) 

where 77^ " — nondiagonal complex 3x3 matrices = 1,2,3). As mentioned in 
the Introduction, in order to define the quark fields mass eigenstates the untary mixing 
matrix V u %^ should be introduced in the Lagrangian terms of the charged Higgs boson 
interaction with quarks 

MteP tfy „£ H + 1 Mu TV u j h - (hQ ) 
r- u L v u i >dJ a R n + a L v i d jU R n . (oyj 

\/2v \/2vtgp 
If we extract the universal phase factor from the mixing matrix elements V u i & — > e 



op 



dJ 



u,a J 



v, 



the Yukawa interaction terms take the form 



^ , L , V u ^ d R H + + u R H~, (60) 

so we can identify the universal phase as the relative phase £ of the SU{2) doublets. 
The structure of this sort, however, does not look like the weak charged current sector 
mixing matrix, where the universal complex factor is not suitable to describe the effects 
of CP violation in meson decays. 
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Appendix 

The decay width h{ — > 77 can be written as 



T(hi -> 77) 



Ml 



a 



\SUM ht )\ 2 + \P?(M hi )\ : 



2567T 3 v 2 

where the scalar and the pseudoscalar factors are given by [Q ESj 



(61) 



S]{M hi ] 



2 E N c Q}g s hif - f — F sf (r if ) 
f=b,t,xt,xf mf 

/j=il,t2,bl,fe2,Tl,T2 /; 



■Pb(Tiif±) 



2 E N cQ)9h i 



f=b,t,xf,xf 



If 



(62) 



Tj X = M|./4m 2 , A^c = 3 for squarks and Nc = 1 for stau and chargino, respectively The 
vertex factors g^fj can De easily extracted from Table 6, where we list also the triple 
vertices with hi and gauge bosons. The threshold corrections induced by the exchanges 
of gluinos and charginos |24| 12*7] are not included in the following calculation. 

The factors F s f, F p f, F F 1 |2E| are expressed by means of the dimensionless function 
f(r) 

Fsf(r) = t- 1 [1 + (1 - t- 1 )^)] , F pf (r) = r" 1 f(r) , (63) 
F (r) = r-M-l+T-VW], F 1 (t) = 2 + 3^ + 3^(2 - r" 1 )/^) , 



with an integral repesentation 
1 fidy 



f(r) 



2 Jo y 



In [l-4ry(l -y)\ 



arcsin 2 ( v / r) : 



in 



T<1, 
T > 1 . 



(64) 



QCD corrections in the large mass limit can be found in [2E] 



V 



71 



Jl 



1 + 



3tt 



(65) 



Chargino contributions depend on the couplings 



s 

hixtXi 
P 

hixtXi 



V n U 12 GS x + V l2 U n GS 2 , 
V n U 12 GP 1 + V 12 U U GP 2 



(66) 



9 h 
9 h 



1X2 X2 



1X2 X2 



V 2 \U 22 GS\ + V 22 U 2 \ GS 2 . 
V 21 U 22 GP 1 + V 22 U 21 GP 2 



(67) 
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for hi we have GS\ = — sin a an + cos a 0,21, GS2 = cos a an + sin a 021, GP\ 
GP 2 = cos/5 a 3 i, and the matrix elements Uy 



1 M| -/i 2 - 2 m 2 ^ cos 2/5 

l/ 12 _l7 21 __yi+ - 



£ B / M 2 2 -/^- 2m*, cos 2/3 

^ = -^ 11 = ^V 1 ^ 



T/ T/ e A Mf - jg+H cos 2/3 

V2l -- Vl2 -^V + w 



T/ T/ 4 / M|-/i 2 + 2 m 2 ,, cos 2(5 

where 

W = \J (Mf + /i 2 + 2 m 2 w y -4(M 2 -fi-m 2 w sin 2/3) 2 , 
£a = sign(M 2 sin/3 + /i cos/5), Eb = sign(M 2 cos/3 + /x sin/5). 
Chargino masses are given by 

m\ x = i| v /(M 2 -/i 2 ) 2 + 2m&,(l + sin2/3) 

- ^(M 2 + /x 2 ) 2 + 277^(1 - sin 2/3)|, 

4 = 1 -^(Ml-^y + 2m 2 w (l + sm2(3) 

+ ^(M 2 + /i 2 ) 2 + 2m 2 w (l - sin 2/3)). 
Sfermion contributions depend on the couplings 

a = (1, 2, a), /3, 7 = L, R, % = (h, h 2 , h 3 ) = (1, 2, 3) j, k — 1,2, 

{ smfl/e**/ cos 9 f ) ' 
r i/V" = _ r <Pif*f sin a + r fc/* / CQS a> 
r 2/*/ = r^-^cosa + F^sina:, 



where 



v^2 V -ih b (spA b + cppi*) 
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Y<t>i~b*b 
pat*t 

Y<Pit*t 
Y<l>2t*t 

par *f 

p0lf *f 
p</>2T*f 



L/i* 4* 



-\h b \ 2 vc() + \ (gl + \gfj vcp 
-j= 2 hA b 

~\ (pi + \gt) vs p ±h$n 

J_( 2/l*( C/3 A* + S/3 /i) 

^2 V -ih t {cpA t + S/3/i*) 

\h t \ 2 vs p + \ (gl - \gf) vsp 
1 

v^2 

_W ih* T (spA* T + cpii) 

y/2 \ -ihrispAr + cpfi*) 

h r \ 2 vcf3 + \ (gl -gDvcp 



\9l) vcp 



J_U* A* 



TsMf -|^t| 2 ^S/3+^ 2 ^ 



L/i* 4* 



In these formulas h t>b>T are real variables. 
Sfermion masses are given by 



m?,r> 

9(0l,2 



where 



m 



qL 



+ m 2 + c 2p m 2 z (T* - Q q s 2 w ), 



m \n = M l 3 + m \ + c 2 pm 2 z Q q s 2 W) 
a q m q = h q v q (A q - /i*R q )/V2, 
m lL = M L + + C2 P m 2 z (s 2 w - 1/2), 



Mi 



C2f3m 2 z s 2 w , 



aitrii = h T v\ {A^ — n* tan (3) / \/2. 

Here the Yukawa couplings of quarks h q , q — t,b, R = U,D, T\ = 
Qb = —1/3, R b = tg/3 = v 2 /v 1 , R t = ctg/5, the mixing angles are 

-\a q{ i)\m q{ i) 



(77) 



-T b z = 1/2, Q t = 2/3, 



cos 9 



9(0 



sin# 



5(0 



v / K-(0^ m ?-(0i )2 + |a ^ |2m ^)' 
2 2 

m q(l)L~ m q (l) 1 
VK-(i)i" m 9-(0l )2+|a ^ |2m 9(')' 
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(79) 



Charged Higgs boson contribution depends on the effective triple self-couplings g , 
which can be written as 

9h+ h- Ax = ~\ — \ — (4(s Q • c/ + c a ■ s p 3 )s 2p m H 2 a 2 i - 8cp 2 sp +a a 21 Reii 2 12 

-8sp +a sp 4 a 21 Re^\ 2 - 8cp +a c (3 2 a 11 Re^\ 2 - 8cp +a sp i a 11 Re^\ 2 

-op 2 's 2f 3 2 'sf3 +a a 2 iReX 6 v 2 + 4c /3 2 s a _ /3 a 2 iReA 6 t; 2 

+4cg 4 s Q ,sg 3 aiiReA6'y 2 + AcaC^s^auReXav 2 
+4:Cp 2 s a s paiiReX^v 2 + 4c a c / 3 5 anReA 6 v 2 

+4(c a • c/ - s a ■ Sf3 3 )s 2/3 m h 2 an - s 2/3 2 sp +a S/3 2 a 21 ReX 7 v 2 

—4:S a ^/3S/3 2 a 2 iRe\ 7 v 2 + Acp 2 s a s p 5 anRe\ 7 v 2 

—4c a C[3Sp 2 aiiRe\ 7 v 2 — 4s Q s / 3 3 anReA 7 f 2 

-Ac a cp 3 sp A a u Re\ 7 v 2 + 4:Cp_ a s 2p 2 m H ± 2 a 21 - 4s 2/3 2 s a _ /3 m H ± 2 ai 1 
-s 2 p 3 Sp +a m A 2 a 21 - cp +a s 2 p 3 m A 2 a u - s 2/3 3 s p +a a 2 iRe\ 5 v 2 
-Cf3 +a s 2f3 3 a u Re\ 5 v 2 + 8cp 3 sp 3 a 31 Im\ 5 v 2 - 8cp 2 sp A a 31 Ira\ 6 v 2 
— 8cp 4 sp 2 a 31 Im\ 7 v 2 ) 



1 X 

h- h 2 = -7, — o — ms a - C/3 3 + c a - s/3 3 )s 2/3 m H 2 a 22 - 8c/3 2 sp +a a 22 Renl 2 

Z S 2/3 Z ■ V 

4 2 2 2 4 2 

—8sp +a Sfs a 22 Re/i 12 — 8cp +a Cp a 12 Refi 12 — 8cp +a Sp ai 2 Refi 12 

-cp 2 s 2p 2 sp +a a 22 Re\ 6 v 2 + 4c /3 2 s a _ /3 a 22 ReA 6 t; 2 

+Acp A s a Sp 3 a 12 Re\ 6 v 2 + Ac a Cp 3 Sp A ai 2 Re\ 6 v 2 
+4cp 2 s a spai 2 Re\ 6 v 2 + 4c a cg 5 a 12 ReA 6 -u 2 

+4(c a • C/3 3 - s a ■ Sp 3 )s 2 pm h 2 a l2 - s 2 p 2 Sp +a Sp 2 a 22 Re\ 7 v 2 

—As a -pSp 2 a 22 Re\ 7 v 2 + Acp 2 s a Sp 5 a 12 Re\ 7 v 2 

— Ac a Cps p 2 a\ 2 Re\ 7 v 2 — 4s Q s / 3 3 ai 2 ReA7f 2 

-4c a cp 3 Sf3 4 a 12 Re\ 7 v 2 + Acp_ a s 2 p 2 m H ± 2 a 22 - 4:S 2(3 2 s a _ /3 m H ± 2 a 12 
-s 2 p 3 sp +a m A 2 a 22 - cp+ a s 2 p 3 m A 2 a 12 - s 2 p 3 sp +a a 22 Re\ 5 v 2 
-cp +a s 2 p 3 ai 2 Re\ 5 v 2 + 8cp 3 sp 3 a 32 Im\ 5 v 2 - 8cp 2 sp A a 32 Im\ e v 2 
-8cp A sp 2 a 32 Im\ 7 v 2 ) 

1 1 

9h+ h- h A = -7, — o — ( 4 ( s « ' c /3 3 + c « • S/3 3 )s 2/3 m H 2 a 23 - 8cp 2 sp +a a 23 Re/i 2 12 

2 S 2 p z ■ V 
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— 8sp +a sp a 23 Refi 12 — 8cp +a Cp ai 3 Re/i 12 — 8cp +Q Sp ai 3 Re/i 12 

-cp 2 s 2 p 2 sp +a a 23 ReX 6 v 2 + Acp 2 s a _pa 23 ReX 6 v 2 

+Acp A s a sp 3 a 13 ReX 6 v 2 + 4c Q ,c /3 3 s /3 4 ai 3 ReA 6 t> 2 
+4cp 2 s a s pdi 3 ReX 6 v 2 + Ac a cp 5 ai 3 ReX 6 v 2 

+4(c Q • c/ - s Q • s /3 3 )s 2/ 3m /l 2 ai 3 - s 2 p 2 Sp +a Sp 2 a 23 ReX 7 v 2 

-As a -pSp 2 a 23 ReX 7 v 2 + 4c« 2 ,s Q ,s / 3 5 ai 3 ReA 7 t> 2 

— Ac a cps p 2 ay 3 ReX 7 v 2 — 4 1 s Q ,Sg 3 ai 3 ReA 7 t> 2 

-4c Q c / g 3 s /3 4 ai 3 ReA 7 t; 2 + 4c /3 __ Q ,s 2/3 2 m // ± 2 a 23 - As 2 p 2 s a _pm H ± 2 a 13 
-S2/3 3 s / 3 +a m A 2 a 23 - c /3+a s 2/ 3 3 m^ 2 ai 3 - s 2 p 3 sp +a a 23 ReX 5 v 2 
-cp +a s 2 p 3 a 13 ReX 5 v 2 + 8c /3 3 s /3 3 a 33 ImA 5 t; 2 - 8c /3 2 s /3 4 a 33 ImA 6 v 2 
-8cp 4: sp 2 a 33 ImX 7 v 2 ) 

This representation uses the mass basis for CP even/odd Higgs fields (h,H,A) then 
rotated by matrix in the three-dimensional (h, H, A) isospace, and for this reason 
includes rrih, ran, itla and m#± of the CP conserving limit, calculated with one-loop 
MSSM corrections from the squark sector. In this sense the vertices above are MSSM 
effective one-loop Higgs self-interaction vertices. If the imaginary parts in these vertices 
are set to zero they are reduced to the self-interaction vertices of the CP conserving limit, 
when m h , ran, vha an d mn± are the masses of physical states. Various extremal cases 
(decoupling limits) are clearly seen. Equivalent representation of the triple couplings can 
be written in the Aj basis (see details on the representations in mass and Aj basis in j5J). 
For example 

3 

q , = -« V a a i q , 

a=l 

where 

9 1H+H - = ReAA 5 spcp c a+ p - ReAA 6 c a s 2 p cp 

+ ReAA 6 s a sp + ReAA 7 Cp (s a Sp cp - c a Up - 2 Spjj 

- ReAA 6 s a s 2 p cp - 2 s a sp cp Ai + 2 c a sp c 2 p A 2 

- cp s a A 3 + c a sp A 3 — c a c 2 p sp A4 + cp s a s 2 p A4, 
9 2H + H ~ = RsAAs spCp s a +p + 2 ReAAg c a sp Cp ReAA 6 c a sp 

- ReAA 6 s a s 2 pCp - ReAX 7 cp (c a spCp + s a (c| — 2s|)j 
+2 c a s 2 p cpXx + 2 s a sp c 2 pX 2 + c a cpX 3 

+ s a spX 3 — c a s 2 pCpX^ — s a spc 2 pXi, 
9 3H+H - = c p ImAA 7 - sp cp ImAA 5 + s 2 ^ ImAA 6 

In this representation the scalar masses of the CP conserving limit do not explicitly 
participate. The magnitude of the coupling g H + H - hl is shown in Figffl 
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Figure 1: Triple Higgs boson interaction vertex g H + hl (GeV) vs the phase Arg(fiA) at 
the parameter values M SU sy =500 GeV, tg/3 =5, A t)6 =1000 GeV, // =2000 GeV. Dashed 
line m H ± =300 GeV, solid line m#± =200 GeV. 



The decay width hi — > (7(7 has the form 



r(Ai -> ##) 



32tt 3 



where 



(80) 



/= fe >* /j=tl, 6,61,62 



= 

and QCD i^-f actors are 



E 9h if f—Fpf(Tif) 
f=b,t " l f 



1 a s (Ml) /95 _ 7 



7T 



4 6 



K 9 A = 1 + 



«s«)/97 7 



7T 



4 6 



52) 



iVV = 5 is the number of quark flavors with masses less than . 

The decay width of Higgs boson to the two fermions hi — > // can be written as 

Ncg}m hl p% { (s a a 2 i - c a a u ) 2 jr + ctg 2 /? a 2 31 , / = u, c, t, 
h ^ff~ I (c a a 21 - s a a u fi + tg 2 ( 3af !1 , f = b,d,s,e, h,t, 



57T 



(83) 



where (3k = 1 — 4/c, k 



9f = 2m^ anc ^ =3 (1) for quarks (leptons) 
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In the following Table 5 we list the Higgs boson masses m^, m/, 2 , rrih z which are 
calculated using the effective \ parameters (fT3 ]) -([T§ J) . Section 2, and the mass term diag- 
onalization method described in Section 3.1. The decay widths T hl _^ gg , T fcl ^ 77 (unprimed) 
include only the leading one-loop contributions of t, b quarks and W ± bosons. For an 
illustration of the sensitivity of , rrih 2 , and their decay widths to the values of Aj 
we computed Higgs boson masses m 'h 1 ,h 2 ,h 3 an d the leading one-loop decay widths Y' hi ^ gg , 
r/ii->77 (i ncmc le t, b and W contributions only) using the effectife potential parametriza- 
tion with both the one- loop and two-loop contributions to Aj from the paper [3| . Finally, 
the decay widths T'l n _ >gg) are found using the effective parameters (fT3*j) -(fT^ j) and 

taking into accout all possible one- loop fermion (t, b), gauge boson W ± , sfermion (t, 
b), chargino and charged Higgs boson contributions, with ^-factors introduced in the 
expressions for decay widths. 

Table 5 contains also the output of the CPsuperH [T7j package and the FeynHiggs [T%] 
package with the input parameter values taken the same as used in our parameter set. 
The two-loop evaluation in the CPsuperH and the one-loop evaluation in the FeynHiggs 
2.1beta has been performed. Note that physical Higgs bosons Hi, H2, H3 of the CPsuperH 
and FeynHiggs are evaluated in the way that is technically different from the construction 
of our mixed states hi, h 2 , h 3 , however a difference of numbers (which is from several 
percent to 40% in the majority of cases) is caused mainly by theoretical uncertainties 
of the effective two-doublet potential representation, not by different definitions of the 
Higgs boson eigenstates in the generic basis of scalar doublets, as demonstrated explicitly 
in section 3.1. 

In Fig.5 and Figs. 6-9 we show the variation of the light Higgs boson mass and the 
variations of T(hi — > gg), T(hi — > 77) decay widths in different regions of the parameter 
space (<p, m H ±, A tjb , /i, tg/3). At the parameter set ( 0, 300 GeV, 1000 GeV, 2000 GeV, 
5 ) the decay widths of hi to 77 and gg are not far from the decay widths of the SM 
Higgs boson with run =120 GeV. Largest sensitivity of the widths to the charged Higgs 
mass is observed. At m#± around 200 GeV (Fig.6a, Fig. 8a) we observe the suppression 
of the branchings of hi to gg and 77 of more than 10 times at (p ~ it, which takes place 
in CPsuperH and FeynHiggs at higher masses of m H ± around 300 GeV. 

Our approach is algorithmized in the form of the model in CompHEP 41.10 format 
|3()j . where the symbolic expressions for vertices are a starting level for calculation of 
the complete tree- level sets of diagrams with the following cross section/decay width 
calculations and the generation of unweighted events. 
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Table 5: Higgs boson mases and their two-particle decay widths. The parameter set ci-EM(x a z) = 0.7812-10 2 , 
a s (m z ) =0.1172, G F = 1.174- 10- 5 GeV" 2 , m b =3 GeV, tg/3 =5, MgvsY = 500 GeV, \A t \ = \A b \ = A = 1000 GeV, 
= 2000 GeV, m H ± = 300 GeV. Our results together with CPsuperH [17] and FeynHiggs [18] with options 2003011100 
(the one- loop regime). , T denote our results with the at one-loop, m' h , V our results with the two-loop terms 
introduced to A^, T/V" are the decay widths in our case when sparticles are not involved/included. 
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Table 6: Vertex factors of hi, h 2 , h 3 . This is a part of the complete set of vertices 
generated by LanHEP package |31. 
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Figure 2: Neutral Higgs boson masses h,H,A versus m#± and the trilinear parameters 
A t , At in the CP conserving limit. Solid line denotes the rrih mass, short dashed - tua, 
long-dashed - m H . (a) tg/3 =5, M SUSY =0.5 TeV, A t = A b = fi =0; (b) tg/J =5, 
M sc/sy =0.5 TeV, A* = A fe =0.9 TeV, // = -1.5 TeV; (c) tg/3 =5, M 5[/5 y =0.5 TeV, 
m^± =220 GeV, // =0, A t = A b ; (d) tg/3 =5, M sc/sy =0.5 TeV, m^± =220 GeV, /i = -2 
TeV, A t = A 6 . 
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Figure 3: (a) Neutral Higgs boson masses, (b),(c),(d) the matrix elements versus the 
phase ip = arg(/iA) at the parameter values tg/? =5, m H ± =180 GeV, M S usy =0.5 TeV, 
A t = A b =1 TeV, /i =2 TeV. Solid line denotes i = 1, long dashed - i = 2 and short 
dashed - i — 3. 
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Figure 4: (a) Neutral Higgs boson masses, (b),(c),(d) the matrix elements versus the 
phase tp = axg(fiA) at the parameter values tg(3 =5, m H ± =300 GeV, M S usy =0.5 TeV, 
A t = A b =1 TeV, /i =2 TeV. Solid line denotes i = 1, long dashed - i = 2 and short 
dashed - i — 3. 
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Figure 5: Light Higgs boson mass m,h 1 (<f) (GeV) vs arg(yUv4) in various regions of 
the MSSM parameter space. Horizontal dotted lines indicate the hi mass in the CP- 
conserving limit (m^j = rrih). (a) tg/3 =5, M S usy =0.5 TeV, A t = A b =1 TeV, /i =2 TeV. 
Solid line m H ± =180 GeV, dashed m H ± =250 GeV. Thin solid line denotes m h ({p). (b) 
tg/3 =5, m H ± =300 GeV, \i =2 TeV, solid line A t = A b =-1.2 TeV, dashed A t = A b =1.3 
TeV. (c) tg/3 =5, m H ± =300 GeV, A t = A b =1 TeV, solid line /i =-1.6 TeV, dashed 
/i =0.7 TeV (d) n =2 TeV, m H ± =300 GeV, A t = A b =1 TeV, solid line tg/3 =5, dashed 
tg/3 =40. 
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Figure 6: The decay width T(hi -> x 10 4 (GeV) at M SU sy =500 GeV. Dotted lines 
show I\ in the CP conserving limit <p = 0. Thin solid or dashed lines are only for SM 
contributions, thick solid or dashed lines show SM and sparticle contrbutios with K-factor 
included, (a) tg/3 =5, A t = A b =1 TeV, fi =2 TeV, solid line m H ± =190 GeV, dashed 
m H± =300 GeV, (b) tg/3 =5, m H ± =300 GeV, // =2 TeV, solid line A t = A b =-1.1 TeV, 
dashed A t = A b =1.1 TeV, (c) tg/3 =5, m^± =300 GeV, A t = A b =1 TeV, solid line 
/i =0.2 TeV, dashed // =1.2 TeV, (d) /i =2 TeV, m H ± =300 GeV, A = A b =1 TeV, solid 
line tg/3 =5, dashed tg/3 =40 
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Figure 7: The decay width T(h x -> x 10 4 (GeV) at M 5C/5 y =500 GeV. Dotted lines 
show Th in the CP conserving limit ip = 0. Thin solid or dashed lines are only for SM 
contributions, thick solid or dashed lines show SM and sparticle contributions with K- 
factor included, (a) tg/? =5, A t = A b = 1 TeV, /i =2 TeV, solid line <p = tt/2, dashed 
<p = it, (b) tg/3 =5, m^± =300 GeV, /i =2 TeV, solid line <p = tt/2, dashed <p = it, (c) 
tg/3 =5, m#± =300 GeV, A t = Ab =1 TeV, solid line ip = tt/2, dashed <p — ir, (d) fi =2 
TeV, m H ± =300 GeV, A t = At, =1 TeV, solid line ip = ir/2, dashed <p = n. 
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Figure 8: The decay width T(h! -> 77) x 10 6 (GeV) at M S ijsy =500 GeV. Dotted 
lines show I\ in the CP conserving limit (p — 0. Thin solid or dashed lines are only for 
SM contributions, thick solid or dashed lines show SM and sparticle contributions with 
J-factor included, (a) tg/3 =5, A t = A b =1 TeV, fi =2 TeV, solid line m H ± =190 GeV, 
dashed m H ± =300 GeV, (b) tg/3 =5, m H ± =300 GeV, /1 =2 TeV, solid line A t = A b =-1.1 
TeV, dashed A t = A b =1.1 TeV, (c) tg/3 =5, m H ± =300 GeV, A t = A b =1 TeV, solid line 
/1 =0.2 TeV, dashed \i =1.2 TeV, (d) /1 =2 TeV, m^± =300 GeV, A t = A b =1 TeV, solid 
line tg/3 =5, dashed tg/3 =40 
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Figure 9: The decay width T{h x -> 77) x 10 6 (GeV) at M SU sy =500 GeV. Dotted 
lines show Th in the CP conserving limit ip = 0. Thin solid or dashed lines denote the 
SM contributions, thick solid or dashed lines show SM and sparticle contributions with 
J-factor included, (a) tg/3 =5, A t = A b =1 TeV, \i =2 TeV, solid line ip = ir/2, dashed 
p = 7r, (b) tg/3 =5, m H ± =300 GeV, \i =2 TeV, solid line ip = n/2, dashed <p = n, (c) 
tgP =5, rriH± =300 GeV, A t = Ai, =1 TeV, solid line <£> = 7r/2, dashed p = n, (d) fi =2 
TeV, =300 GeV, A t = Ab =1 TeV, solid line <p = ir/2, dashed <p — ir. 
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